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ABSTRACT 

In a warm dark matter (WDM) cosmology, the first objects to form at z > 20 are one dimensional 
filaments with mean length on the order of the WDM free-streaming scale. Gao and Theuns recently 
claimed by using high-resolution hydrodynamic simulations that the eventual collapse of these WDM 
filaments along their longest axes may seed the supermassive black holes that power high-z quasars. In 
this picture, it is supposed that the high-z quasar luminosity function should reflect how abundant the 
WDM filaments are in the early universe. We derive analytically the mass function of early-universe 
filaments with the help of the Zel'dovich approximation. Then, we determine the rate of its decrease in 
the mass section corresponding to the free streaming scale of a WDM particle of mass m v . Adjusting 
the value of m v , we fit the slope of the analytic model to that of the high-z quasar luminosity function 
measured from the Sloan Digital Sky Survey DR3. A new WDM constraint from this feasibility study 
is found to be consistent with the lightest super-symmetric partner. 
Subject headings: cosmology: theory — large-scale structure of universe 



1. INTRODUCTION 

The large-scale features of the observed universe are 
strikingly consistent with the theoretical predictions 
based on the cold dark matter model. The combined 
analyses of the recent data from the observations of cos- 
mic microwave background (CMB) , galaxy power spec- 
trum and Type la supernovae (e.g.. lDunklev et aLll2009l . 
and references therein) have been capable of measuring 
the key cosmological parameters that characterize the 
CDM model with surprisingly high precision. This has 
opened an era of precision cosmology, echoing the tri- 
umph of the CDM model. 

Nevertheless, the status of the CDM model as the stan- 
dard paradigm has been shaking currently in both obser- 
vational and theoretical perspectives. Observations have 
reported several mismatches between the predictions of 
the CDM model and the real phenomena on galactic 
and subgalactic scales. For instance, the abundance 
of galactic satellites, the slope of the inner core of the 
dark halo density profiles, and the degree of void empti- 
ness have exhibi ted apparent conflicts between theory 
and observat ion (jKlvpin et all 1 19991: iMoore et alJll99a 
lPeeblesl[200lT ). Very recently. iDisnev et al.l (|2008f) mea- 
sured the cross-correlations between different galaxy's 
observables and showed that only a single parameter suf- 
fices to explain the complex structures of galaxies. This 
phenomenon is inconsistent with the CDM picture where 
the galaxy formation is driven by hierarchical merging 
process. Although it is still inconclusive whether the re- 
ported mismatches really indicate the failure of the CDM 
model or they are just by-products of the selection effects 
and/or limitations of the observational techniques, it has 
certainly led to an emergence of alternative models. 

A theoretical challenge comes from the fact that the ex- 
tension of the standard particle physics favors the warm 
dark matter (WDM) rather than CDM. The WDM dif- 
fers from CDM in a respect that the WDM particles 
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have non- negligible free streaming scale, l v . On scales 
larger than l v , WDM behave just like CDM. On the other 
hand, on scales less than l u , WDM do not cluster but 
free stream due to their large velocity dispersions unlike 
CDM which cluster on all scales. It is this property that 
enables the WDM model to overcome the observational 
mismatches faced by the CDM model. The list of realis- 
tic WDM candidates that have so far proposed includes 
gravitin os dBertone et al.ll2005l and references therein), 
axions ( Ba ttve fc Shellardl 11994ft , and sterile neutrinos 
(|Dodelson fc Widrowi 11994 ). each of which has distinct 
characteristic mass range. By constraining the WDM 
particle mass, it may be possible to rule out a certain 
WDM candidate. In fa ct there have been already such 
attempts: For instance, ISeljak et al.l (|2006l ) constrained 
the particle mass of WDM by using the Lya forest power 
spectrum and sug gested tha t the sterile neutrinos be 
ruled out (see also lViel et aHl2006() 

In a WDM universe, the objects that will first con- 
dense out in the initial density inhomogeneities are not 
zero dimensional clum ps but one dimensional filaments. 
iGao fc Theunsl (|2007l ) recently studied the formation and 
properties of WDM filaments with high-resolution hy- 
drodynamic simulations, and proposed that the early- 
universe filaments would seed the supermassive black 
holes that power the quasars at z > 6 through their 
eventual collapse along the longest axes. If this sce- 
nario is true, then, the following two issues should be 
related to each other: How luminous and abundant the 
first generation high-z quasars are; and how massive and 
abundant the early-universe WDM filaments are. In 
statistical terms, the luminosity function of the high-z 
quasars should be related to the mass function of the 
early-universe filaments. Given that the mass function 
of the early-universe filaments depends on the particle 
mass of WDM (see §2), it implies the possibility of using 
the high-z quasar luminosity function as a new WDM 
constraint. 

The idealistic way to study the mass function of early- 
universe filaments and its dependence on the WDM par- 
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ticle mass is to use the hydrodynamics simulations. How- 
ever, given the computational costs as well as the current 
resolution-limit of hydrodynamic simulations, it is ade- 
quate and quite necessary to consider an analytic ap- 
proach as the first guideline. Here, we evaluate ana- 
lytically the abundance of early-universe filaments and 
constrain the WDM particle mass by fitting its slope in 
the WDM particle free streaming scale to that of the 
high- z quasar luminosity function (jFan et al.ll2001l ) de- 
termined from the Sloan Digital S ky Survey Data Re- 
lease 3 (SDSS) (jTrump et alJl2006Q . For the key cosmo- 
logical parameters other than the WDM particle mass , 
we assume a WMAP5 cosmology (jDunklev et al.ll2009t ) 
throughout this paper. 

The plan of this paper is as follows. In §2, an analytic 
model for the abundance of early universe filaments is 
presented. In §3, a constraint on the particle mass of 
WDM is provided by fitting the analytic model to the 
observational data from SDSS. In §4, the achievements 
as well as the caveats of our work are discussed. 

2. THE ABUNDANCE OF EARLY UNIVERSE FILAMENTS 

The Zel'dovich approximation predicts generically and 
in the mo st simple way the formation of early-universe 
filaments (IZel'Dovichl |1970h . According to this model, 
the mass density p of a given region can be expressed in 
terms of the Lagrangian quantities as 



P = 



[1 - £>(z)Ai][l - D(z)X 2 ][l - D(z)X 3 ] ' 



(1) 



where p is the mean background density, Ai, A2, A3 (in a 
decreasing order) are the three eigenvalues of the linear 
deformation tensor, and D(z) is a linear growth factor. 
Equation (JTJ) predicts that p will diverge at Ai = 1/D(z). 
That is, the first collapse will occur when the largest 
eigenvalue reaches a threshold value, A c (z) = 1/D(z). 
The sign of the other two eigenvalues, A2 and A3, at the 
moment of the first collapse will determine the dimension 
of a collapsed object. A one-dimensional filament will 
form if A2 > and A3 < 0, while the formation of a two- 
dimensional sheet will occur for the case A2 < and A3 < 
0. Since we are interested in the formation of "filaments" 
that are shown t o be most abundant in the early universe 
filled by WDM (|Gao fc Theund 120071) . we consider only 
the filament condition, \i = 1/D(z), A 2 > and A 3 < 0. 

The fractional volume F(M, z) occupied by those La- 
grangian regions which will condense out early-universe 
filaments at z on mass scale M can be expressed as 

F(M, z) = p[Ai > A c (z), A 2 > 0, A 3 < : a M ], (2) 

where A c (z) = 1/D(z) and gm is the rms density fluc- 
tuation related to the dimensionless linear WDM power 
spectrum A^(fc) as 



(T 



.2 _ 
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Al(k)W 2 (k,M)d\nk, 



(3) 



where W(k,M) 
mass scale M. 



represents a sharp k-space filter on 
We adopt the approximation formula 



o fBode et all ((200l . A^(fc) = A 2 (yfc)|T^| where A 2 (fc) 
is the CDM linear power spectrum. The functional form 
of the WDM transfer function, T% is given as 
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where 

a = 0.048(^/0.4) - 15 (V0.65) 1 - 3 (keV/™,) 1 - 15 (1.5/.g l/ ) - 29 , 

7 = 1.2, m v is the WDM particle mass and g v is the num- 
ber of degrees of freedom that equals 3/4 for fermions. 

The fractional volume in equation Q can be ob- 
tained by integrating the joint probability density dis- 
tri bution of the three e igenvalues that was first derived 
by iDoroshkevichl (|1970D as 
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where 1\ = Ai + A2 + A3 and I2 — A1A2 + A2A3 + 
A3A1. In the spirit of the Press-Schechter theory 
(|Press &: Schechtenll974 hereafter PS), we evaluate the 
number density of early-universe filaments collapsed at 
mass scale M as 



n F (M, z) 



P_ 

M 



dM 



F(M, z) 



(7) 



(4) 



where the redshift for the formation of the early- 
uni verse filaments is set at z = 23, in accordance with 
iGao fc Theunsl (|2007[ ). In the original PS theory which 
dealt with the mass function of the clumps, equation l[7]l 
had to be multiplied by an additional normalization con- 
stant to account fo r the occurrence of the cloud-in-clouds 
(|Bond et alj|1991l) . But, for our case the cloud-in-clouds 
will not occur since the WDM particles just free-stream 
in the early-universe filaments forming no substructures. 
Henceforth, we do not have to account for the cloud- 
in-cloud problem in deriving the mass function of early 
universe filaments. 

The cumulative mass function n F (> M) of early uni- 
verse filaments can be obtained by integrating equation 
as n F (> M) = n F (M)dM. Since our focus is 
not on the amplitude of the mass function but mainly 
on its slope, we rescale n F (> M) to have unity as M 
goes zero. Figure Q] plots the rescaled cumulative mass 
function np(> M) at z — 23 versus M/Mq for the case 
of vn v — IkeV (open circles) and m v = lOkeV (open 
squares). Here M$ corresponds to the free streaming 
mass scale of a given WDM particle. As can be seen, 
7ip(> M) begins to decrease as M becomes larger than 
Mq. The rate of its decrease at M > M depends sen- 
sitively on the value of the WDM particle mass m v : It 
decreases more rapidly with M for the case of smaller 
value of m„. It is because the less massive WDM par- 
ticles have smaller free-stream scales. To quantify the 
rate of the decrease of rip(> M) at M > Mo, we fit 
n F {> M) to a power-law M^ 1 (i.e., n F {M) cx M 13 ) in 
the mass range where 0.1 < n F (> M)/n F (> Mo) < 1. 
In the higher mass section where n F {> F) drops expo- 
nentially with M, we expect bas ically no early univers e 
filaments. As demonstrated by IGao &; Theunsl (|2007h . 
only those early-universe filaments with mass compara- 
ble in order of magnitude to the WDM free-streaming 
scale can form which are in turn related to the first- 
generation supermassive black holes that power the ob- 
served high-z luminous quasars. The higher mass sec- 
tion where n F (> M)/n F (> Mq) < 0.1 correspond to the 
scales order-of-magnitude larger than the free streaming 
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Fig. 1. — The rescaled cumulative mass function of early-universe 
filaments for the two different cases of WDM particle mass (m„ = 1 
and fOkeV as open circles and open squares, respectively). The 
value Mo corresponds to the free-streaming mass scale of a given 
WDM particle. The solid lines represent the power-law fitting in 
the mass range where O.f < np(> M)/np(> Mq) < 1 for the two 
cases. 

scale of a given WDM particles. The solid lines in Fig. Q] 
represent the power-law fitting of the cumulative mass 
function of early universe filaments. As can be seen, the 
slope (3 indeed depends sensitively on the value of m v . 

3. RELATION TO HIGH-Z QUASAR LUMINOSITY 
FUNCTION 

The discovery of a supermassive black hole with mass 
3 x IO^Mq in qua sar SDSS J1148+5251 at z = 6.41 
(|Willott et al. 2003) has raised a crucial question of what 
seeded such a supermassive black hole at that early epoch 
(|Podsi adlowsk i ~et al.ll2003 | . and re f erence s therein). As 
introduced in §1, iGao fc Theund (|2007l ) have studied 
the formation of early-universe filaments using high- 
resolution hydrodynamic simulations for a WDM cosmol- 
ogy and suggested that the ultimate collapse of the WDM 
filaments along the longest axes should induce vigorous 
collisions between gas clouds and stars due to their high 
densities, which would in turn seed the formation of su- 
permassive black holes that power quasars at z > 6. In 
subsequent evolution, the black holes would grow active 
by accreting baryonic gases and dark matter particles 
that constituted their parent filaments. In this picture, 
the number density of active black holes at high redshifts 
would be related to that of early universe filaments. 

The luminosity functio n <fi (L) o f high -z quasars at 
z > 4 was determined bv iFan et al.l (|2001l ). who showed 
that (j) (L) can be well approximated as a power law with 
high-end slope in the range of — 2.5±0.25. If these high-z 
quasars are indeed powered by the first-generation super- 
massive black holes formed through the ultimate collapse 
of the WDM filaments, the high-end slope of the high-z 
quasar luminosity function should be consistent with the 
power-law slope of the mass function of early-universe 
filaments determined in §2. Varying the value of m„, we 
determine the power-law slopes of the mass function of 
early-universe filaments (3, which are shown in Figure O 
The lcr and 2a range of the high-end slope of the high-z 
quasar luminosity function cj) (L) measured from SDSS 



Fig. 2. — The power-law slope of the mass function of early uni- 
verse filaments versus WDM particle mass. The horizontal dashed 
and dotted lines correspond to the la and 2tr ranges of the high- 
end slope o f high-z quasar luminosity function determined from 
SDSS DR3 ijFan et aLII2u0ir i. respectively. 

DR3 are shown as horizontal dashed and dotted lines, re- 
spectively. The comparison indicates that the two slopes 
agree with each other when m„ is a few keV. 

4. DISCUSSION AND CONCLUSION 

Motivated by the recent heuristic work of 

iCxao fc Theunl (pOOl . we have constructed an an- 
alytic model for the mass function of the early-universe 
filaments in a WDM universe. Then we test the 
possibility of using it as a new WDM constraint. As a 
feasibility study, we compare the high-end slope of the 
mass function of the early-universe filaments to that of 
the high-z quasar luminosity function determined from 
SDSS DR3 and find a new WDM constraint, m v ~ a 
few keV. This preliminary re sult looks consisten t with 
the super-symmetric picture (jBertone et al.l [2005). The 
advantage of the mass function of the early-universe 
filament as a WDM constraint lies in the fact that the 
model is purely analytical, free of any fitting or nuisance 
parameter. It depends only on the WDM particle 
mass other than the key cosmological parameters. In 
addition, when the analytic model is compared to the 
observational result, it is not required to account for 
any complicated baryon physics since we consider only 
the slope of the mass function. 

Yet, there are a couple of simplified assumptions on 
which our model is based. First, the redshifts of the 
quasars should be higher than z > 6 for a more fair 
comparison. In this work, we just used the redshift 
range z > 4, assuming that the quasars at z > 4 cor- 
respond to the first generation. To be more realistic, it 
will be required to determine the luminosity function of 
the quasars at z > 6, but there are only small number 
of quasars have so far been observed at th at high red- 
shift. For example, Vestcrg arrd et all (|2008l ) determined 
the supermassive black hole mass function at z > 4.7 us- 
ing a tens of quasars and found the high-end slope in a 
much larger range of — 1.9 ± 1.7, which obviously suffers 
from large uncertainty due to the small number statistics. 
Anyway, we would like to mention clearly that the true 
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high-end slope of the high-z quasar luminosity function 
is currently significantly more uncertain than the data 
used here. 

Second, the analytic model for the mass function of 
early universe filaments has to be refined. Our model is 
the simplest approximation based on the Zel'dovich ap- 
proximation and the Press-Schechter theory. The valid- 
ity of our model has to be tested numerically before using 
it in practice. Third, the evolution of the early-universe 
filaments has not been taken into account properly for 
the evaluation of the mass function. In reality, before 
their eventual collapse at z > 6, the early-universe fila- 
ments may undergo fragmentation due to the tidal effect 
from the surrounding matter. To account for the evolu- 



tionary effects, it will be anyway required to refine the 
model with the help of the hydrodynamic simulations. 
Our future work is in this direction. 

As a final conclusion, our model and the preliminary 
result from our feasibility study has provided a proof of 
concept that the high-z quasar luminosity function will 
be in principle a useful WDM probe. 
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